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1. Introduction 


ABSTRACT 


We compare the performance of a recently proposed empirical climate model based on astronomical 
harmonics against all CMIP3 available general circulation climate models (GCM) used by the IPCC 
(2007) to interpret the 20th century global surface temperature. The proposed astronomical empirical 
climate model assumes that the climate is resonating with, or synchronized to a set of natural 
harmonics that, in previous works (Scafetta, 2010b, 2011b), have been associated to the solar system 
planetary motion, which is mostly determined by Jupiter and Saturn. We show that the GCMs fail to 
reproduce the major decadal and multidecadal oscillations found in the global surface temperature 
record from 1850 to 2011. On the contrary, the proposed harmonic model (which herein uses cycles 
with 9.1, 10-10.5, 20-21, 60-62 year periods) is found to well reconstruct the observed climate 
oscillations from 1850 to 2011, and it is shown to be able to forecast the climate oscillations from 1950 
to 2011 using the data covering the period 1850-1950, and vice versa. The 9.1-year cycle is shown to be 
likely related to a decadal Soli/Lunar tidal oscillation, while the 10-10.5, 20-21 and 60-62 year cycles 
are synchronous to solar and heliospheric planetary oscillations. We show that the IPCC GCM’s claim 
that all warming observed from 1970 to 2000 has been anthropogenically induced is erroneous because 
of the GCM failure in reconstructing the quasi 20-year and 60-year climatic cycles. Finally, we show 
how the presence of these large natural cycles can be used to correct the IPCC projected anthropogenic 
warming trend for the 21st century. By combining this corrected trend with the natural cycles, we show 
that the temperature may not significantly increase during the next 30 years mostly because of the 
negative phase of the 60-year cycle. If multisecular natural cycles (which according to some authors 
have significantly contributed to the observed 1700-2010 warming and may contribute to an 
additional natural cooling by 2100) are ignored, the same IPCC projected anthropogenic emissions 
would imply a global warming by about 0.3-1.2 °C by 2100, contrary to the IPCC 1.0-3.6 °C projected 
warming. The results of this paper reinforce previous claims that the relevant physical mechanisms that 
explain the detected climatic cycles are still missing in the current GCMs and that climate variations at 
the multidecadal scales are astronomically induced and, in first approximation, can be forecast. 

© 2011 Elsevier Ltd. All rights reserved. 


The astronomical harmonic model assumes that the climate 
system is resonating with or is synchronized to a set of natural 


Herein, we test the performance of a recently proposed 
astronomical-based empirical harmonic climate model (Scafetta, 
2010b, in press) against all general circulation climate models 
(GCMs) adopted by the IPCC (2007) to interpret climate change 
during the last century. A large supplement file with all GCM 
simulations herein studied plus additional information is added 
to this manuscript. A reader is invited to look at the figures 
depicting the single GCM runs there reported to have a feeling 
about the performance of these models. 
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frequencies of the solar system. The synchronicity between solar 
system oscillations and climate cycles has been extensively 
discussed and argued in Scafetta (2010a,b, 2011b), and in the 
numerous references cited in those papers. We used the velocity 
of the Sun relative to the barycenter of the solar system and a 
record of historical mid-latitude aurora events. It was observed 
that there is a good synchrony of frequency and phase between 
multiple astronomical cycles with periods between 5 and 100 
years and equivalent cycles found in the climate system. We refer 
to those works for details and statistical tests. The major 
hypothesized mechanism is that the planets, in particular Jupiter 
and Saturn, induce solar or heliospheric oscillations that induce 
equivalent oscillations in the electromagnetic properties of the 
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upper atmosphere. The latter induces similar cycles in the cloud 
cover and in the terrestrial albedo forcing the climate to oscillate 
in the same way. The soli/lunar tidal cyclical dynamics also 
appears to play an important role in climate change at specific 
frequencies. 

This work focuses only on the major decadal and multidecadal 
oscillations of the climate system, as observed in the global 
surface temperature data since AD 1850. A more detailed discus- 
sion about the interpretation of the secular climate warming 
trending since AD 1600 can be found in Scafetta and West (2007) 
and in Scafetta (2009) and in numerous other references there 
cited. About the millenarian cycle since the Middle Age a discus- 
sion is present in Scafetta (2010a) where the relative contribution 
of solar, volcano and anthropogenic forcing is also addressed, and 
in the numerous references cited in the above three papers. Also 
correlation studies between the secular trend of the temperature 
and the geomagnetic aa-index, the sunspot number and the solar 
cycle length address the above issue and are quite numerous: for 
example Hoyt and Schatten (1997), Sonnemann (1998), and Thejll 
and Lassen (2000). Thus, a reader interested in better under- 
standing the secular climate trending topic is invited to read 
those papers. In particular, about the 0.8 °C warming trending 
observed since 1900 numerous empirical studies based on the 
comparison between the past climate secular and multisecular 
patterns and equivalent solar activity patterns have concluded 
that at least 50-70% of the observed 20th century warming could 
be associated to the increase of solar activity observed since the 
Maunder minimum of the 17th century: for example see Scafetta 
and West (2007), Scafetta (2009), Loehle and Scafetta (2011), 
Soon (2009), Soon et al. (2011), Kirkby (2007), Hoyt and Schatten 
(1997), Le Mouél et al. (2008), Thejll and Lassen (2000), Weihong 
and Bo (2010), and Eichler et al. (2009). Moreover, Humlum et al. 
(2011) noted that the natural multisecular/milennial climate 
cycles observed during the late Holocene climate change clearly 
suggest that the secular 20th century warming could be mostly 
due to these longer natural cycles, which are also expected to cool 
the climate during the 21th century. A similar conclusion has 
been reached by another study focusing on the multisecular and 
millennial cycles observed in the temperature in the central- 
eastern Tibetan Plateau during the past 2485 years (Liu et al., 
2011). For the benefit of the reader, in Section 7 in the supple- 
ment file the results reported in two of the above papers are very 
briefly presented to graphically support the above claims. 

It is important to note that the above empirical results contrast 
greatly with the GCM estimates adopted by the IPCC claiming that 
more than 90% of the warming observed since 1900 has been 
anthropogenically induced (compare Figures 9.5a and b in the IPCC 
report which are reproduced in Section 4 in the supplement file). In 
the above papers it has been often argued that the current GCMs 
miss important climate mechanisms such as, for example, a modula- 
tion of the cloud system via a solar induced modulation of the cosmic 
ray incoming flux, which would greatly amplify the climate sensi- 
tivity to solar changes by modulating the terrestrial albedo (Scafetta, 
2011b:; Kirkby, 2007; Svensmark, 1998, 2007; Shaviv, 2008). 

In addition to a well-known decadal climate cycle commonly 
associated to the Schawbe solar cycle by numerous authors (Hoyt 
and Schatten, 1997), several studies have emphasized that the 
climate system is characterized by a quasi bi-decadal (from 
18 year to 22 year) oscillation and by a quasi 60-year oscillation 
(Stockton et al., 1983; Currie, 1984; Cook et al., 1997; Agnihotri 
and Dutta, 2003; Klyashtorin et al., 2009; Sinha et al., 2005; 
Yadava and Ramesh, 2007; Jevrejeva et al., 2008; Knudsen et al., 
2011; Davis and Bohling, 2001; Scafetta, 2010b; Weihong and Bo, 
2010; Mazzarella and Scafetta, 2011; Scafetta, in press). For 
example, quasi 20-year and 60-year large cycles are clearly 
detected in all global surface temperature instrumental records 


of both hemispheres since 1850 as well as in numerous astro- 
nomical records. There is a phase synchronization between these 
terrestrial and astronomical cycles. As argued in Scafetta (2010b), 
the observed quasi bidecadal climate cycle may also be around a 
21-year periodicity because of the presence of the 22-year solar 
Hale magnetic cycle, and there may also be an additional 
influence of the 18.6-year soli/lunar nodal cycle. However, for 
the purpose of the present paper, we can ignore these corrections 
which may require other cycles at 18.6 and 22 years. In the same 
way, we ignore other possible slight cycle corrections due to the 
interference/resonance with other planetary tidal cycles and 
with the 11-year and 22-year solar cycles, which are left to 
another study. 

About the 60-year cycle it is easy to observe that the global 
surface temperature experienced major maxima in 1880-1881, 
1940-1941 and 2000-2001. These periods occurred during the 
Jupiter/Saturn great conjunctions when the two planets were 
quite close to the Sun and the Earth. This events occur every three 
J/S synodic cycles. Other local temperature maxima occurred 
during the other J/S conjunctions, which occur every about 20 
years: see Figures 10 and 11 in Scafetta (2010b), where this 
correspondence is shown in details through multiple filtering of 
the data. Moreover, the tides produced by Jupiter and Saturn in 
the heliosphere and in the Sun have a period of about 
0.5/(1/11.86—1/29.45) +10 years plus the 11.86-year Jupiter 
orbital tidal cycles. The two tides beat generating an additional 
cycle at about 1/(2/19.86—1/11.86)=61 years (Scafetta, in 
press). Indeed, a quasi 60-year climatic oscillations have likely 
an astronomical origin because the same cycles are found in 
numerous secular and millennial aurora and other solar related 
records (Charvatova et al., 1988; Komitov, 2009; Ogurtsov et al., 
2002; Patterson et al., 2004; Yu et al., 1983; Scafetta, 2010a,b, 
2011b; Mazzarella and Scafetta, 2011). 

A 60-year cycle is even referenced in ancient Sanskrit texts 
among the observed monsoon rainfall cycles (Iyengar, 2009), a fact 
confirmed by modern monsoon studies (Agnihotri and Dutta, 2003). 
It is also observed in the sea level rise since 1700 (Jevrejeva et al., 
2008) and in numerous ocean and terrestrial records for centuries 
(Klyashtorin et al., 2009). A natural 60-year climatic cycle associated 
to planetary astronomical cycles may also explain the origin of 60- 
year cyclical calendars adopted in traditional Chinese, Tamil and 
Tibetan civilizations (Aslaksen, 1999). Indeed, all major ancient 
civilizations knew about the 20-year and 60-year astronomical 
cycles associated to Jupiter and Saturn (Temple, 1998). 

In general, power spectrum evaluations have shown that 
frequency peaks with periods of about 9.1, 10-10.5, 20-22 and 
60-63 years are the most significant ones and are common 
between astronomical and climatic records (Scafetta, 2010b, in 
press). Evidently, if climate is described by a set of harmonics, it 
can be in first approximation reconstructed and forecast by using 
a planetary harmonic constituent analysis methodology similar to 
the one that was first proposed by Lord Kelvin (Thomson, 1881; 
Scafetta, in press) to accurately reconstruct and predict tidal 
dynamics. The harmonic constituent model is just a superposition 
of several harmonic terms of the type 


N 
F(t) =Ao+ X Ai costœit+ Q), @ 


i=1 


whose frequencies œ; are deduced from the astronomical theories 
and the amplitude A; and phase ¢, of each harmonic constituent 
are empirically determined using regression on the available data, 
and then the model is used to make forecasts. Several harmonics 
are required: for example, most locations in the United States use 
computerized forms of Kelvin’s tide-predicting machine with 
35-40 harmonic constituents for predicting local tidal amplitudes 
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(Ehret, 2008), so a reader should not be alarmed if many harmonic 
constituents may be needed to accurately reconstruct the climate 
system. 

Herein we show that a similar harmonic empirical methodology 
can, in first approximation, reconstruct and forecast global climate 
changes at least on a decadal and multidecadal scales, and that this 
methodology works much better than the current GCMs adopted by 
the IPCC (2007). In fact, we will show that the IPCC GCMs fail to 
reproduce the observed climatic oscillations at multiple temporal 
scales. Thus, the computer climate models adopted by the IPCC 
(2007) are found to be missing the important physical mechanisms 
responsible for the major observed climatic oscillations. An impor- 
tant consequence of this finding is that these GCMs have seriously 
misinterpreted the reality by significantly overestimating the 
anthropogenic contribution, as also other authors have recently 
claimed (Douglass et al., 2007; Lindzen and Choi, 2011; Spencer and 
Braswell, 2011). Consequently, the IPCC projections for the 21st 
century should not be trusted. 


2. The IPCC GCMs do not reproduce the global surface 
temperature decadal and multidecadal cycles 


Fig. 1 depicts the monthly global surface temperature anomaly 
(from the base period 1961-1990) of the Climatic Research Unit 
(HadCRUT3) (Brohan et al., 2006) from 1850 to 2011 against an 
advanced general circulation model average simulation (Hansen 
et al., 2007), which has been slightly shifted downward for visual 
convenience. The chosen units are the degree Celsius in agree- 
ment with the climate change literature referring to temperature 
anomalies. The GISS ModelE is one of the major GCMs adopted by 
the IPCC (2007). Here we study all available climate model 
simulations for the 20th century collected by Program for Climate 
Model Diagnosis and Intercomparison (PCMDI) mostly during the 
years 2005 and 2006, and this archived data constitutes phase 3 
of the Coupled Model Intercomparison Project (CMIP3). These 
GCMs, use the observed radiative forcings (simulations 
“tas:20c3m”) adopted by the IPCC (2007). All GCM simulations 
are depicted and analyzed in Section 2 of the supplement file 
added to this paper. These GCM simulations cover a period that 


1 T 


may begin during the second half of the 19th century and end 
during the 21th century. The following calculations are based on 
the maximum overlapping period between each model simula- 
tion and the 1850-2011 temperature period. The CMIP3 GCM 
simulations analyzed here can be downloaded from Climate 
Explorer web-site: see the supplement file for details. 

A simple visual inspection suggests that the temperature 
presents a quasi 60-year cyclical modulation oscillating around 
an upward trend (Scafetta, 2010b, 2011b; Loehle and Scafetta, 
2011). In fact, we have the following 30-year trending patterns: 
1850-1880, warming; 1880-1910, cooling; 1910-1940, warming; 
1940-1970, cooling; 1970-2000, warming; and it is almost 
steady or presents a slight cooling since 2001 (2001-2011.5 
rate= — 0.46 + 0.3 °C/ century). Other global temperature recon- 
structions, such as the GISSTEM (Hansen et al., 2007) and the GHCN- 
Mv3 by NOAA, present similar patterns (see Section 1 in the 
supplement file). Note that GISSTEM/1200 presents a slight warm- 
ing since 2001 (2001-2011.5 rate=+047 +0.3°C/ century), 
which appears to be due to the GISS poorer temperature sampling 
during the last decade for the Antarctic and Arctic regions that were 
artificially filled with a questionable 1200 km smoothing methodol- 
ogy (Tisdale, 2010). However, when a 250 km smooth methodology 
is applied, as in GISSTEM/250, the record shows a slight cooling 
during the same period (2001-2011.5 rate=—0.16 +0.3°C/ 
century). HadCRUT data has much better coverage of the Arctic 
and Southern Oceans that GISSTEM and, therefore, it is likely more 
accurate. Note that CRU has recently produced an update of their 
SST ocean record, HadSST3 (Kennedy et al., 2011), but it stops in 
2006 and was not merged yet with the land record. This new 
corrected record presents an even clearer 60-year modulation than 
the HadSST2 record because in it the slight cooling from 1940 to 
1970 is clearer (Mazzarella and Scafetta, 2011). 

Indeed, the 60-year cyclicity with peaks in 1940 and 2000 
appears quite more clearly in numerous regional surface tem- 
perature reconstructions that show a smaller secular warming 
trending. For example, in the United States (D’Aleo, 2011), in the 
Arctic region (Soon, 2009), in several single stations in Europe and 
other places (Le Mouél et al., 2008) and in China (Soon et al., 
2011). In any case, a 60-year cyclical modulation is present for 
both the Norther and Southern Hemisphere and for both Land and 
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Fig. 1. Global surface temperature (top, http://www.cru.uea.ac.uk/cru/data/temperature/) and GISS ModelE average simulation (bottom). The records are fit with Eq. (5). 
Note also the large volcano eruption signatures that appear clearly overestimated in the GCM’s simulation. 
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Ocean regions (Scafetta, 2010b) even if it may be partially hidden 
by the upward warming trending. The 60-year modulation 
appears well correlated to a recently proposed solar activity 
reconstruction (Loehle and Scafetta, 2011). 

The 60-year cyclical modulation of the temperature from 1850 
to 2011 is further shown in Fig. 2 where the autocorrelation 
functions of the global surface temperature and of the GISS 
ModelE average simulation are compared. The autocorrelation 
function is defined as 


N-t T T 
m DaD O tay 


Jenura] 


where T is the average of the N-data long temperature record and 
t is the time-lag. The autocorrelation function of the global 
surface temperature (Fig. 2A) and of the same record detrended 
of its quadratic trend (Fig. 2B) reveals the presence of a clear 
cyclical pattern with minima at about 30-year lag and 90-year lag, 
and maxima at about 0-year lag and 60-year lag. This pattern 
indicates the presence of a quasi 60-year cyclical modulation in 
the record. Moreover, because both figures show the same pattern 
it is demonstrated that the quadratic trend does not artificially 
creates the 60-year cyclicity. On the contrary, the GISS ModelE 
average simulation produces a very different autocorrelation 
pattern lacking any cyclical modulation. Fig. 2C shows the auto- 
correlation function of the two records detrended also of their 
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Fig. 2. Autocorrelation function (Eq. (2)) of the global surface temperature and of the GISS ModelE average simulation: [A] original data; [B] data detrended of their 
quadratic fit; [C] the 60-year modulation is further detrended. Note the 60-year cyclical modulation of the autocorrelation of the temperature with minima at 30-year and 
90-year lags and maxima at 0-year and 60-year lags, which is not reproduced by the GCM simulation. Moreover, the computer simulation presents an autocorrelation peak 
at 80-year lag related to a pattern produced by volcano eruptions, which is absent in the temperature. See Section 5 in the supplement file for further evidences about the 
GISS ModelE serious overestimation of the volcano signal in the global surface temperature record. 
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60-year cyclical fit, and the climatic record appears to be 
characterized by a quasi 20-year smaller cycle, as deduced by 
the small but visible quasi regular 20-year waves, at least up to a 
time-lag of 70 years after which other faster oscillations with a 
decadal scale dominate the pattern. On the contrary, the auto- 
correlation function of the GCM misses both the decadal and bi- 
decadal oscillations and again shows a strong 80-year lag peak, 
absent in the temperature. The latter peak is due to the quasi 
80-year lag between the two computer large volcano eruption 
signatures of Krakatoa (1883) and Agung (1963-1964), and to the 
quasi 80-year lag between the volcano signatures of Santa Maria 
(1902) and El Chichón (1982). Because this 80-year lag autocor- 
relation peak is not evident in the autocorrelation function of the 
global temperature we can conjecture that the GISS ModelE is 
significantly overestimating the volcano signature, in addition to 
not reproducing the natural decadal and multidecadal tempera- 
ture cycles: this claim is further supported in Section 5 of the 
supplement file. 
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A similar qualitative conclusion applies also to all other GCMs 
used by the IPCC, as shown in Section 2 of the supplement file. 
The single GCM runs as well as their average reconstructions 
appear quite different from each other: some of them are quite 
flat until 1970, others are simply monotonically increasing. 
Volcano signals often appear overestimated. Finally, although 
these GCM simulations present some kind of red-noise variability 
supposed to simulate the multi-annual, decadal and multidecadal 
natural variability, a simple visual comparison among the simula- 
tions and the temperature record gives a clear impression that the 
simulated variability has nothing to do with the observed tem- 
perature dynamics. In conclusion, a simple visual analysis of the 
records suggests that the temperature is characterized 10-year, 
20-year and 60-year oscillations that are simply not reproduced 
by the GCMs. This is also implicitly indicated by the very smooth 
and monotonically increasing pattern of their average reconstruc- 
tion depicted in the IPCC figure SPM.5 (see Section 4 in the 
supplement file). 
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Fig. 3. [A] Maximum Entropy Method (MEM) with M=N/2 (solid) poles and the Lomb Periodogram (dash) of the HadCRUT3 global surface temperature monthly sampled 
from 1850 to 2011 (see Section 3 of the supplement file for details and explanations). The two techniques produce the same peaks, but MEM produces much sharper peaks. 
The major four peaks are highlighted in the figure. [B] As above for the HadCRUT, GISSTEM/250 and GHCN-Mv3 global surface temperature records during the period 
1880-2011: see Section 1 in the supplement file. Note that the spectra are quite similar, but for GISSTEM the cycles are somehow slightly smoother and smallers than for 


the other two sequences, as the bottom curves show. 
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Fig. 3A and B shows two power spectra estimates of the 
temperature records based on the Maximum Entropy Method 
(MEM) and the Lomb periodogram (Press et al., 2007). Four major 
peaks are found at periods of about 9.1, 10-10.5, 20-21 and 
60-62 years: other common peaks are found but not discussed 
here. Both techniques produce the same spectra. To verify 
whether the detected major cycles are physically relevant and 
not produced by some unspecified noise or by the specific 
sequences, mathematical algorithms and physical assumptions 
used to produce the HadCRUT record, we have compared the 
same double power spectrum analysis applied to the three 
available global surface temperature records (HadCRUT3, GIS- 
STEM/250 and GHCN-Mv3) during their common overlapping 
time period (1880-2011): see also Section 1 in the supplement 
file. As shown in the figures the temperature sequences present 
almost identical power spectra with major common peaks at 
about 9.1, 10-10.5, 20 and 60 years. Note that in Scafetta (2010b), 
the relevant frequency peaks of the temperature were determined 
by comparing the power spectra of HadCRUT temperature records 
referring to different regions of the Earth such as those referring 
to the Northern and Southern hemispheres, and to the Land and 
the Ocean. So, independent major global surface temperature 
records present the same major periodicities: a fact that further 
argues for the physical global character of the detected 
spectral peaks. 

Note that a methodology based on a spectral comparison of 
independent records is likely more physically appropriate than 
using purely statistical methodologies based on Monte Carlo 
randomization of the data, that may likely interfere with weak 
dynamical cycles. Note also that a major advantage of MEM is that 
it produces much sharper peaks that allow a more detailed 
analysis of the low-frequency band of the spectrum. Section 5 in 
the supplement file contains a detailed explanation about the 
number of poles needed to let MEM to resolve the very-low 
frequency range of the spectrum: see also Courtillot et al. (1977). 

Because the temperature record presents major frequency 
peaks at about 20-year and 60-year periodicities plus an appar- 
ently accelerating upward trend, it is legitimate to extract these 
multidecadal patterns by fitting the temperature record (monthly 
sampled) from 1850 to 2011 with the 20 and 60-year cycles plus a 
quadratic polynomial trend. Thus, we use a function f(t)+p(t) 
where the harmonic component is given by 


27(t—T,) i 27(t—T2) 
fom=Q cos| 60 | -Ca cos 0 |. (3) 
and the upward quadratic trending is given by 
p(t) = P2*(t—1850)* + Py(t—1850)+ Po. (4) 


The regression values for the harmonic component are 
Cı =0.10 + 0.01 °C and C2 = 0.040 + 0.005 °C, and the two dates 
are Tı = 2000.8 + 0.5 AD and T3 = 2000.8 + 0.5 AD. For the quad- 
ratic component we find Pp =—0.30+0.2 °C, Pı = —0.0035 + 
0.0005 °C/year and P, = 0.000049 + 0.000002 °C/year?. Note that 
the two cosine phases are free parameters and the regression 
model gives the same phases for both harmonics, which suggests 
that they are related. Indeed, this common phase date approxi- 
mately coincides with the closest (to the sun) conjunction 
between Jupiter and Saturn, which occurred (relative to the 
Sun) on June/23/2000 ( ~ 2000.5), as better shown in Scafetta 
(2010b). 

It is important to stress that the above quadratic function p(t) 
is just a convenient geometrical representation of the observed 
warming accelerating trend during the last 160 years, not outside 
the fitting interval. Another possible choice, which uses two linear 
approximations during the periods 1850-1950 and 1950-2011, 
has also be proposed (Loehle and Scafetta, 2011). However, our 


quadratic fitting trending cannot be used for forecasting purpose, 
and it is not a component of the astronomical harmonic model. 
Section 4 will address the forecast problem in details. 

It is possible to test how well the IPCC GCM simulations 
reproduce the 20 and 60-year temperature cycles plus the 
upward trend from 1850 to 2011 by fitting their simulations 
with the following equation: 


m(t) = ax0.10 cos pe 


2n(t—2000.8) 
20 


+ bx0.040 cos| | +cxp(t)+d, (5) 
where a, b, c and d are regression coefficients. Values of a, b and c 
statistically compatible with the number 1 indicate that the 
model well reproduces the observed temperature 20 and 60-year 
cycles, and the observed upward temperature trend from 1850 to 
2011. On the contrary, values of a, b and c statistically incompa- 
tible with 1 indicate that the model does not reproduce the 
observed temperature patterns. 

The regression values for all GCM simulations are reported in 
Table 1. Fig. 4 shows the values of the regression coefficients a, b 
and c for the 26 climate model ensemble-mean records and all fail 
to well reconstruct both the 20 and the 60-year oscillations found 
in the climate record. In fact, the values of the regression 
coefficients a and b are always well below the optimum value 
of 1, and for some model these values are even negative. The 
average among the 26 models is a=0.30+0.22 and 
b=0.35 + 0.42, which are statistically different from 1. This result 
would not change if all available single GCM runs are analyzed 
separately, as extensively shown in Section 2 of the supplement 
file. 

About the capability of the GCMs of reproducing the upward 
temperature trend from 1850 to 2011, which is estimated by the 
regression coefficient c, we find a wide range of results. The 
average is c=1.11 + 0.50, which is centered close to the optimum 
value 1. This result explains why the multi-model global surface 
average simulation depicted in the IPCC figures 9.5 and SPM.5 
apparently reproduces the 0.8 °C warming observed since 1900. 
However, the results about the regression coefficient c vary 
greatly from model to model: a fact that indicates that these 
GCMs usually also fail to properly reproduce the observed upward 
warming trend from 1850 to 2011. 

Table 1 and the tables in Section 2 in the supplement file also 
report the estimated reduced y? values between the measured 
GCM coefficients dm, bm and Cm (index “m” for model) and the 
values of the same coefficients ar, br and cr (index “T”’ for 
temperature) estimated for the temperature. The reduced x? 
(chi square) values for three degree of freedom (that is three 
independent variables) are calculated as 
(bm—br)* 


Ab? + Aby. 


(Cm—Cr ? 
Ac2 + Act 


2_1 (am-an? 
~ 3 | Aa? + Aas 


x (6) 
where the A values indicate the measured regression errors. We 
found x? > 1 for all models: a fact that proves that all GCMs fail to 
simultaneously reproduce the 20-year, 60-year and the upward 
trend observed in the temperature with a probability higher than 
99.9%. This y? measure based on the multidecadal patterns is 
quite important because climate changes on a multidecadal scale 
are usually properly referred to as climate changes, and a climate 
model should at least get these temperature variations right to 
have any practical economical medium-range planning utility 
such as street construction planning, agricultural and industrial 
location planning, prioritization of scientific energy production 
research versus large scale applications of current very expensive 
green energy technologies, etc. 
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Table 1 

Values of the regression parameters of Eq. (5) obtained by fitting the 25 IPCC (2007) climate GCM ensemble-mean estimates. #1 refers to the ensemble average of the GISS 
ModelE depicted in Fig. 1a; #2-#26 refers to the 25 IPCC GCMs. Pictures and analysis concerning all 95 records including each single GCM run are shown in Section 2 in the 
supplement file that accompanies this paper. The optimum value of these regression parameters should be a=b=c=1 as presented in the first raw that refers to the 
regression coefficients of the same model used to fit the temperature record. The last column refers to a reduced 7? test based on three coefficients a, b and c: see Eq. (6). 
This determines the statistical compatibility of the regression coefficient measured for the GCM models and those observed in the temperature. It is always measured a 
reduced x? > 1 for three degrees of freedom, which indicates that the statistical compatibility of the GCMs with the observed 60-year, 20-year temperature cycles plus the 
secular trending is less than 0.1%. These GCM regression values are depicted in Fig. 4: the regression coefficients for each available GCM simulation are reported in the 
supplement file. The y? test in the first line refers to the compatibility of the proposed model in Eq. (3) relative to the ideal case of a=b=c=1 that gives a reduced 
X? =0.21 which imply that the statistical compatibility of Eq. (3) with the temperature cycles plus the secular trending is about 90%. The fit has been implemented using 
the nonlinear least-squares (NLLS) Marquardt-Levenberg algorithm. 


# Model a (60-year) b (20-year) c (trend) d (bias) 7? (abc) 
Temp 1.03 + 0.05 0.99 + 0.12 1.01 + 0.02 0.00 + 0.01 0.21 
1 GISS ModelE 0.25 + 0.03 0.90 + 0.08 0.80 + 0.01 0.08 + 0.01 89 
2 BCC CM1 0.63 + 0.03 0.69 + 0.09 0.54 + 0.02 0.08 + 0.01 109 
3 BCCR BCM2.0 0.29 + 0.05 0.06 + 0.11 0.40 + 0.02 0.08 + 0.01 202 
4 CGCM3.1 (T47) 0.35 + 0.03 —0.28 + 0.07 2.02 + 0.01 0.40 + 0.01 753 
5 CGCM3.1 (T63) 0.11 + 0.05 0.05 + 0.11 2.07 + 0.02 0.40 + 0.01 536 
6 CNRM CM3 —0.01 + 0.07 —0.27 + 0.18 2.02 + 0.03 0.39 + 0.01 322 
7 CSIRO MK3.0 0.30 + 0.04 —0.12 + 0.11 0.48 + 0.02 0.08 + 0.01 176 
8 CSIRO MK3.5 —0.19 + 0.04 —0.19 + 0.10 1.38 + 0.02 0.25 + 0.01 197 
9 GFDL CM2.0 0.44 + 0.05 0.90 + 0.12 1.12 + 0.02 0.21 + 0.01 28 
10 GFDL CM2.1 0.37 + 0.07 0.75 + 0.17 1.37 + 0.03 0.26 + 0.01 53 
11 GISS AOM 0.22 + 0.03 —0.14 + 0.06 1.10 + 0.01 0.22 + 0.01 93 
12 GISS EH 0.48 + 0.04 0.96 + 0.11 0.80 + 0.02 0.14 + 0.01 43 
13 GISS ER 0.47 + 0.04 0.80 + 0.08 0.90 + 0.02 0.11 + 0.01 31 
14 FGOALS g1.0 0.10 + 0.09 —0.15 + 0.21 0.28 + 0.03 0.06 + 0.01 171 
15 INVG ECHAM4 —0.12 + 0.05 0.37 + 0.12 1.34 + 0.02 0.24 + 0.01 138 
16 INM CM3.0 0.30 + 0.07 0.47 + 0.18 1.34 + 0.03 0.24 + 0.01 54 
17 IPSL CM4 0.13 + 0.06 0.05 + 0.14 1.37 + 0.02 0.26 + 0.01 107 
18 MIROC3.2 Hires 0.35 + 0.05 0.92 + 0.12 1.43 + 0.02 0.19 + 0.01 104 
19 MIROC3.2 Medres 0.34 + 0.03 0.76 + 0.09 0.72 + 0.01 0.14 + 0.01 104 
20 ECHO G 0.58 + 0.04 0.16 + 0.10 0.98 + 0.02 0.18 + 0.01 26 
21 ECHAM5/MPI-OM 0.19 + 0.04 0.31 + 0.09 0.70 + 0.02 —0.02 + 0.01 104 
22 MRI CGCM 2.3.2 0.31 + 0.03 0.03 + 0.07 1.36+0.01 0.27 + 0.01 149 
23 CCSM3.0 0.34 + 0.04 0.43 + 0.10 1.29 + 0.02 0.24 + 0.01 76 
24 PCM 0.77 + 0.05 0.49 + 0.12 1.00 + 0.02 0.16 + 0.01 7 
25 UKMO HADCM3 0.28 + 0.05 0.56 + 0.11 0.94 + 0.02 0.18 + 0.01 42 
26 UKMO HADGEM1 0.52 + 0.04 0.63 + 0.10 1.05 + 0.02 0.20 + 0.01 24 
Average 0.30 + 0.22 0.35 + 0.41 1.11 + 0.47 0.19 + 0.11 143.8 


2.5 ~ra aAA 

: : : : : coeff. a for the 60-year cycle —e— 
coeff. b for the 20-year cycle -+-*—4 
; ; : : ; coeff. c for the upward trend -—-*— 
2L : wi zig! : i a=b=c=1 for the upward trend — 


coeff. a, b, c 


climate model number 


Fig. 4. Values of the regression coefficients a, b and c relative to the amplitude of the 60 and 20-year cycles, and the upward trend obtained by regression fit of the 26 GCM 
simulations of the 20th century used by the IPCC. See Table 1 and Section 2 in the supplement file for details. The result shows that all GCMs significantly fail in 
reproducing the 20-year and 60-year cycle amplitudes observed in the temperature record by an average factor of 3. 
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It is also possible to include in the discussion the two detected 
decadal cycles as 


27(t—T3) 27(t—T4) 
10.44 | Ca cos] 9.07 |. 


g(t) = C3 cos| O) 
A detailed discussion about the choice of the two above periods 
and their physical meaning is better addressed in Section 4. 
Fitting the temperature for the period 1850-2011 gives 
C3 = 0.03 + 0.01 °C, T3 =2002.7 + 0.5 AD, C4= 0.05 + 0.01 °C, 
T4 = 1997.7 + 0.3 AD. It is possible to test how well the IPCC 
GCMs reconstruct these two decadal cycles by fitting their 
simulations with the following equation: 


2n(t—2002.7) 


1044 l +l%0.05 cos| 


n(t) = m(t)+s0.03 cos | r | ; 


9.07 
(8) 


where s and l are the regression coefficients. Values of s and | 
statistically compatible with the number 1 indicate that the 
model well reproduces the two observed decadal temperature 
cycles, respectively. On the contrary, values of s and | statistically 
incompatible with 1 indicate that the model does not reproduce 
the observed temperature cycles. The results referring the average 
model run, as defined above, are reported in Table 2, where it is 
evident that the GCMs fail to reproduce these two decadal cycles 
as well. The average values among the 26 models is s=0.06 + 0.40 
and l=0.34 + 0.37, which are statistically different from 1. In 
many cases the regression coefficients are even negative. The 


Table 2 

Values of the regression parameters s and | of Eq. (8) obtained by fitting the 26 
IPCC (2007) climate GCM ensemble-mean estimates. The fit has been implemen- 
ted using the nonlinear least-squares (NLLS) Marquardt-Levenberg algorithm. 
Note that the two regression coefficients are quite different from the optimum 
values s =!=1, as found for the temperature. The column referring to the reduced 
z% test is based on all five regression coefficients (a, b, c, s and I) by extending 
Eq. (6). Again it is always observed a xy? >1, which indicates incompatibility 
between the GCM and the temperature patterns. The last column indicates the 
RMS residual values between the 4-year average smooth curves of each GCM 
simulation and the 4-year average smooth curve of the temperature: the value 
associated to the first raw (temperature) RMS=0.051 °C) refers to the RMS of the 
astronomical harmonic model that suggests that the latter is statistically 2-5 times 
more accurate than the GCM simulations in reconstructing the temperature record. 


# Model s(10.44-year) 1 (9.1-year) 72 (abcsl) RMS (°C) 
0 Temperature 1.06 + 0.16 0.99 + 0.10 0.15 0.051 
1 GISS ModelE 0.30 + 0.11 0.40 + 0.07 61 0.107 
2 BCC CM1 0.53 + 0.11 0.49 + 0.07 70 0.105 
3 BCCR BCM2.0 —0.11 +0.15 0.06 + 0.09 137 0.158 
4 CGCM3.1 (T47) —0.47 +0.09 0.06 + 0.06 479 0.212 
5 CGCM3.1 (T63) 0.39+0.15 —0.11 +0.09 337 0.220 
6 CNRM CM3 0.22 + 0.24 —0.07 +0.14 202 0.229 
7 CSIRO MK3.0 —0.54+0.14 -—0.01+0.09 128 0.169 
8 CSIRO MK3.5 —0.53 +0.13 0.44 + 0.08 134 0.156 
9 GFDL CM2.0 —0.26 +0.16 0.62 +0.10 25 0.113 

10 GFDL CM2.1 0.13 + 0.23 0.98 + 0.14 34 0.170 

11 GISS AOM 0.19 + 0.09 0.10 + 0.05 73 0.101 

12 GISS EH 0.27 + 0.14 0.66 + 0.09 30 0.106 

13 GISS ER 0.29 + 0.11 0.48 + 0.07 25 0.094 

14 FGOALS g1.0 —0.69 +0.29 0.23 +0.17 111 0.252 

15 INVG ECHAM4 —0.35+0.16 —0.23+0.10 105 0.132 

16 INM CM3.0 —0.15 +0.24 1.01 +0.14 36 0.150 

17 IPSL CM4 0.49 + 0.19 0.48 + 0.11 68 0.137 

18 MIROC3.2 Hires 0.17 + 0.16 0.43 + 0.09 69 0.122 

19 MIROC3.2 Medres 0.24 +0.11 0.47 + 0.07 69 0.106 

20 ECHO G 0.52 + 0.13 0.54 + 0.08 20 0.097 

21 ECHAM5/MPI-OM 0.15 +0.12 —0.09 +0.07 82 0.126 

22 MRI CGCM 2.3.2 0.04+0.10 0.25 + 0.06 103 0.114 

23 CCSM3.0 0.12 + 0.13 0.91 + 0.08 50 0.110 

24 PCM 1.01 + 0.16 0.70 + 0.09 5 0.093 

25 UKMO HADCM3 0.07 +0.15 —0.34 +0.09 49 0.123 

26 UKMO HADGEM1 -—0.46+0.14 0.32 + 0.08 30 0.107 

Average 0.06 + 0.40 0.34 + 0.37 97.39 0.139 


table also includes the reduced xy? (chi square) values for five 
degree of freedom by extending Eq. (6) to include the other two 
decadal cycles. Again, we found 7? > 1 for all models. 

Finally, we can estimate how well the astronomical model 
made of the sum of the four harmonics plus the quadratic trend 
(that is f(t)}+g(t)+p(t)) reconstructs the 1850-2011 temperature 
record relative to the GCM simulations. For this purpose we 
evaluate the root mean square (RMS) residual values between 
the 4-year average smooth curves of each GCM average simula- 
tion and the 4-year average smooth of the temperature curve, and 
we do the same between the astronomical model and the 4-year 
average smooth temperature curve. We use a 4-year average 
smooth because the model is not supposed to reconstruct the fast 
sub-decadal fluctuations. The RMS residual values are reported in 
Table 2. The RMS residual value relative to the harmonic model is 
0.051 °C, while for the GCMs we get RMS residual values from 2 to 
5 times larger. This result further indicates that the geometrical 
model is significantly more accurate than the GCMs in recon- 
structing the global surface temperature from 1850 to 2011. 

The above finding reinforces the conclusion of Scafetta (2010b) 
that the IPCC (2007) GCMs do not reproduce the observed major 
decadal and multidecadal dynamical patterns observed in the 
global surface temperature record. This conclusion does not 
change if the single GCM runs are studied. 


3. Reconstruction of the global surface temperature 
oscillations: 1880-2011 


A regression model may always produce results in a reason- 
able agreement within the same time interval used for its 
calibration. Thus, showing that an empirical model can recon- 
struct the same data used for determining its free regression 
parameters would be not surprising, in general. However, if the 
same model is shown to be capable of forecasting the patterns of 
the data outside the temporal interval used for its statistical 
calibration, then the model likely has a physical meaning. In fact, 
in the later case the regression model would be using construc- 
tors that are not simply independent generic mathematical 
functions, but are functions that capture the dynamics of the 
system under study. Only a mathematical model that is shown to 
be able to both reconstruct and forecast (or predict) the observa- 
tions is physically relevant according the scientific method. 

The climate reconstruction efficiency of an empirical climate 
model based on a set of astronomical cycles with the periods 
herein analyzed has been tested and verified in Scafetta (2010b, in 
press) and Loehle and Scafetta (2011). Herein, we simply sum- 
marize some results for the benefit of the reader and for introdu- 
cing the following section. 

In Figures 10 and 11 in Scafetta (2010b) it is shown that the 
20-year and 60-year oscillations of the speed of the Sun relative to 
the barycenter of the solar system are in a very good phase 
synchronization with the correspondent 20 and 60-year climate 
oscillations. Moreover, detailed spectra analysis has revealed that 
the climate system shares numerous other frequencies with the 
astronomical record. 

In Figures 3 and 5 in Loehle and Scafetta (2011) it is shown 
that an harmonic model based on 20-year and 60-year cycles and 
free phases calibrated on the global surface temperature data for 
the period 1850-1950 is able to properly reconstruct the 20-year 
and 60-year modulation of the temperature observed since 1950. 
This includes a small peak around 1960, the cooling from 1940 to 
1970, the warming from 1970 to 2000 and a slight stable/cooling 
trending since 2000. It was also found a quasi linear residual with 
a warming trending of about 0.66+0.16 °C/century that was 
interpreted as due to a net anthropogenic warming trending. 
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In Scafetta (2011b), it was found that the historical mid- 
latitude aurora record, mostly from central and southern Europe, 
presents the same major decadal and multidecadal oscillations of 
the astronomical records and of the global surface temperature 
herein studied. It has been shown that a harmonic model with 
aurora/astronomical cycles with periods of 9.1, 10.5, 20, 30 and 60 
years calibrated during the period 1850-1950 is able to carefully 
reconstruct the decadal and multidecadal oscillations of the 
temperature record since 1950. Moreover, the same harmonic 
model calibrated during the period 1950-2010 is able to carefully 
reconstruct the decadal and multidecadal oscillations of the 
temperature record from 1850 to 1950. The argument about the 
1850-1950-fit versus 1950-2010-fit is crucial for showing the 
forecasting capability of the proposed harmonic model. This 
property is what distinguishes a mere curve fitting exercise from 
a valid empirical dynamical model of a physical system. This is a 
major requirement of the scientific method. A preliminary phy- 
sical model based on a forcing of the cloud system has been 
proposed to explain the synchrony between the climate system 
and the astronomical oscillations (Scafetta, 2011b). 

The above results have supported the thesis that climate is 
forced by astronomical oscillations and can be partially recon- 
structed and forecasted by using the same cycles, but for an 
efficient forecast there is the need of additional information. This 
is done in the next section. 


4. Corrected anthropogenic projected warming trending and 
forecast of the global surface temperature: period 2000-2100 


Even assuming that the detected decadal and multidecadal 
cycles will continue in the future, to properly forecast climate 
variation for the next decades, additional information is neces- 
sary: (1) the amplitudes and the phases of possible multisecular 
and millennial cycles; (2) the net anthropogenic contribution to 
the climate warming according to realistic emission scenarios. 

The first issue is left to another paper because it requires a 
detailed study of the paleoclimatic temperature proxy reconstruc- 
tions which are relatively different from each other. These cycles 
are those responsible for the cooling periods during the Maunder 
and Dalton solar minima as well as for the Medieval Warm Period 
and the Little Ice Age. So, we leave out these cycles here. 
Considering that we may be at the very top of these longer cycles, 
ignoring their contribution may be reasonable only if our forecast 
is limited to the first decades of the 21st century. However, a 
rough preliminary estimate would suggest that these longer 
cycles may contribute globally to an additional cooling of about 
0.1 °C by 2100 because the millenarian cycle presents an approx- 
imate min-max amplitude of about 0.5-0.7 °C (Ljungqvist, 2010) 
and the top of these longer cycles would occur somewhere during 
the 21st century (Humlum et al., 2011; Liu et al., 2011). Secular 
and millennial longer natural cycles could have contributed about 
0.2-0.3 °C warming from 1850 to 2010 (Scafetta and West, 2007; 
Eichler et al., 2009: Scafetta, 2009, 2010a). 

The second issue is herein explicitly addressed by using an 
appropriate argument that adopts the same GHG emission scenar- 
ios utilized by the IPCC, but correct their climatic effect. In fact, the 
combination of the 20-year and 60-year cycles, as evaluated in Eq. 
(3), should have contributed for about 0.3 °C of the 0.5 °C warming 
observed from 1970 to 2000. During this period the IPCC (2007) 
have claimed, by using the GCMs studied herein, that the natural 
forcing (solar plus volcano) would have caused a cooling up to 0.1- 
0.2 °C (see Figure 9.5 in the IPCC report, which is herein reproduced 
with added comments in Figure S3A in Section 4 in the supplement 
file). As it is evident in the IPCC Figure 9.5b (also shown in the 
supplement file), the IPCC GCM results imply that from 1970 to 


2000 the net anthropogenic forcing contributed a net warming of 
the observed 0.5 °C plus, at most, another 0.2 °C, which had to 
offset the alleged natural volcano cooling of up to — 0.2 °C. A 0.7 °C 
anthropogenic warming trend in this 30-year period corresponds 
to an average anthropogenic warming rate of about 2.3 °C/century 
since 1970. This value is a realistic estimate of the average GCM 
performance because the average GCM projected anthropogenic 
net warming rate is 2.3+0.6°C/century from 2000 to 2050 
according to several GHG emission scenarios (see Figure SPM.5 in 
the IPCC report, which is herein reproduced with added comments 
in Figure S4B in the supplement file). 

On the contrary, if about 0.3 °C of the warming observed from 
1970 to 2000 has been naturally induced by the 60-year natural 
modulation during its warming phase, at least 43-50% of the 
alleged 0.6-0.7 °C anthropogenic warming has been naturally 
induced, and the 2.3°C/century net anthropogenic trending 
should be reduced at least to 1.3 °C/century. 

However, the GCM alleged 0.1-0.2 °C cooling from 1970 to 
2000 induced by volcano activity may be a gross overestimation 
of the reality. In fact, as revealed in Fig. 2, the GCM climate 
simulation presents a strong volcano signature peak at 80-year 
time lag that is totally absent in the temperature record, even 
after filtering. This would imply that the volcano signature in the 
global surface temperature record should be quite smaller and 
shorter than what the GCMs estimate, as empirical studies have 
shown (Lockwood, 2008; Thompson et al., 2009). Section 5 of the 
supplement file shows that the GISS ModelE appears to greatly 
overestimate the long-time signature associated to volcano activ- 
ity against the same signature as estimated by empirical studies. 

Moreover, the observed 0.5 °C warming from 1970 to 2000, 
which the IPCC models associate to anthropogenic GHG plus 
aerosol emissions and to other anthropogenic effects, may also 
be partially due to poorly corrected urban heat island (UHI) and 
land use changes (LUC) effects, as argued in detailed statistical 
studies (McKitrick and Michaels, 2007; McKitrick, 2010). As 
extensively discussed in those papers, it may be reasonable that 
the ~ 0.5 °C warming reported since 1950-1970 in the available 
temperature records has been overestimated up to 0.1-0.2 °C 
because of poorly corrected UHI and LUC effects. Indeed, the land 
warming since 1980 has been almost twice the ocean warming, 
which may be not fully explained by the different heat capacity 
between land and ocean. Moreover, during the last decades the 
agencies that provide the global surface temperature records have 
changed several times the methodologies adopted to attempt to 
correct UHI and LUC spurious warming effects and, over time, 
have produced quite different records (D’Aleo, 2011). Curiously, 
the earlier reconstructions show a smaller global warming and a 
more evident 60-year cyclical modulation from 1940 to 2000 than 
the most recent ones. 

Finally, there may be an additional natural warming due to 
multisecular and millennial cycles as explained in Introduction. In 
fact, the solar activity increased during the last four centuries 
(Scafetta, 2009), and the observed global surface warming during 
the 20th century is very likely also part of a natural and persistent 
recovery from the Little Ice Age of AD 1300-1900 (Scafetta and 
West, 2007; Scafetta, 2009; Loehle and Scafetta, 2011; Soon, 
2009; Soon et al., 2011; Kirkby, 2007; Hoyt and Schatten, 1997; 
Le Mouél et al., 2008; Thejll and Lassen, 2000; Weihong and Bo, 
2010; Eichler et al., 2009; Humlum et al., 2011; Liu et al., 2011): 
see also Section 7 in the supplement file. 

Thus, the above estimated 1.30 °C/century anthropogenic warm- 
ing trending is likely an upper limit estimate. As a lower limit we 
can reasonably assume the 0.66 + 0.16 °C/century, as estimated in 
Loehle and Scafetta (2011), which would be compatible with the 
claim that only 0.2 °C warming (instead of 0.7 °C) of the observed 
0.5 °C warming since 1970 could be anthropogenically induced. 
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This result would be consistent with the fact that according empirical 
studies (Lockwood, 2008; Thompson et al., 2009) the cooling long- 
range effects of the volcano eruptions almost vanished in 2000 (see 
Section 5 in the supplement file) and that the secular natural trend 
could still be increasing. So, from 2000 to 2050 we claim that the 
same IPCC (2007) anthropogenic emission projections could only 
induce a warming trend approximately described by the curve 


q(t) = (0.009 + 0.004)(t—2000). (9) 


There are also two major quasi decadal oscillations with periods 
of about 9.1 year and 10-10.5 year: see Fig. 3. The 9.1-year cycle 
may be due to a Soli/Lunar tidal cycle (Scafetta, 2010b, in press). In 
fact, the lunar apsidal line rotation period is 8.85 years while the 
Soli/Lunar nodal cycle period is 18.6 years. Note that there are two 
nodes and the configuration Sun—Moon-Earth and Sun-Earth- 
Moon are equivalent for the tides: thus, the resulting tidal cycles 
should have a period of about 18.6/2=9.3 year. The two cycles at 
8.85-year and 9.3-year should beat, and produce a fast cycle with 
an average period of 2/(1/8.85+ 1/9.93)=9.07 year that could be 
modulated by a slow cycle with period of 2/(1/8.85—1/9.93)= 
182.9 year. There may also be an additional influence of the half 
Saros eclipse cycle that is about 9 years and 5.5 days. In conclusion, 
the quasi 9.1-year cycle appears to be related to a Soli/Lunar tidal 
cycle dynamics. The 10-10.5-year cycle has been interpreted as 
related to an average cycle between the 0.5/(1/11.862 — 1/29.457)= 
9.93 year Jupiter/Saturn half-synodic tidal cycle and the 11-year 
solar cycle (we would have a beat cycle with period of 2/(1/9.93+ 1/ 
11)=10.44 year). Moreover, a quasi 9.91-year and 10.52-year cycles 
have been found in the natural gravitational resonances of the solar 
system (Bucha et al., 1985; Grandpiere, 1996; Scafetta, 2011b). 

It is possible to include these two cycles in the harmonic 
model using the additional harmonic function Eq. (7) and our final 
model based on 4-frequency harmonics plus two independent 
trending functions is made as 


‘ p(t) if 1850 <t < 2000, 
(=FO+8O+ S p2000)+4(t)_ if 2000 < t < 2100, 


To test the forecasting capability of the g(t) harmonics, the 
f()+g(0)+p() model is calibrated in two complementary periods. 
Note that g(t) is sufficiently orthogonal to f(t)+p(t), so we keep 
f(})+p(t) unchanged for not adding too many free regression 
parameters. Fitting the period 1850-1950 gives C3=0.03+ 
0.01 °C, T3=2003+0.5AD, Cy=0.054+0.01 °C, Tz =1997.5+ 
0.3 AD. Fitting the period 1950-2011 gives C3 = 0.04 + 0.01 °C, 
T = 2002.1 + 0.5 AD, Cy=0.05+0.01 °C, T,= 1998.1 + 0.3 AD. 
Fitting the period 1850-2011 gives C3 = 0.03 + 0.01 °C, T3= 
2002.7 + 0.5 AD, C4 = 0.05 + 0.01 °C, T4 = 1997.7 + 0.3 AD. If the 
decadal period 10.44 year is substituted with a 10 year period for 
1850-2011, we get C= 0.02 + 0.01 °C, T3=2000.4+0.5 AD, 
C4 = 0.04 + 0.01 °C, T4 = 1997.7 + 0.3 AD. 

We observe that all correspondent amplitudes and phases 
coincide within the error of measure, which implies that the 
model has forecasting capability. Moreover, the phase related to 
the 9.1-year cycle presents a maximum around 1997-1998. We 
observe that this period is in good phase with the Soli/Lunar nodal 
dates at the equinoxes, when the Soli/Lunar spring tidal maxima 
are located in proximity of the equator, and the extremes in the 
tidal variance occurs (Sidorenkov, 2005). In fact, each year there 
are usually two solar eclipses and two lunar eclipses, but the 
month changes every year and the cycle repeats every about 
9 years with the moon occupying the opposite node. Thus, 
eclipses occur, within a two week interval, close to the equinoxes 
(around March 20/21 and September 22/23) every almost 9 years. 
Section 6 in the supplement file reports the dates of the solar and 
lunar eclipses occurred from 1988 to 2010 and compares these 


(10) 


dates with the detected 9-year temperature cycle. Two lunar 
eclipses occurred on 24/March/1997 and 16/September/1997, the 
latter eclipse also occurred at the lunar perigee (that is, when the 
Moon is in its closest position to the Earth) so that the line of 
the lunar apsides too was oriented along the Earth-Sun direction 
(so that the two cycles could interfere constructively). Two solar 
eclipses took place almost 9-years later at almost the same dates, 
22/September/2006 (at the lunar apogee) and 19/March/2007 (at 
the lunar perigee). This date matching suggests that the 9.1-year 
cycle is likely related to a Soli/Lunar tidal cycle. Indeed, this cycle 
is quite visible in the ocean oscillations (Scafetta, 2010b) and 
ocean indexes such as the Atlantic Multidecadal Oscillation 
(AMO) and the Pacific Decadal Oscillation (PDO). 

The timing of the 10-10.5-year cycle maximum (2000-2003), 
corresponds relatively well with the total solar irradiance max- 
imum in 2002 (Scafetta and Willson, 2009) and the Jupiter/Saturn 
conjunction around 2000.5 (so that the two cycles could interfere 
constructively). This suggests that this decadal cycle has a solar/ 
astronomical origin. 

The above information is combined in Fig. 5A and B that 
depict: the monthly sampled global surface temperature since 
1850; a 4-year moving average estimates of the same; the 
proposed model given in Eq. (10) with two and four cycles, 
respectively. Finally, for comparison, we plot the IPCC projected 
warming using the average GCM projection estimates, which is 
given by a linear trending warming of 2.3 + 0.6 °C/century from 
2000 to 2050 while since 2050 the projections spread a little bit 
more according to alternative emission scenarios (see figure S4B 
in Section 4 in the supplement file). The two figures are com- 
plementary by highlighting both a low resolution forecast that 
extends to 2100, which can be more directly compared with the 
IPCC projections, and a higher resolution forecast for the next 
decades that may be more important for an immediate econom- 
ical planning, as explained above. 

Fig. 5 clearly shows the good performance of the proposed 
model (Eq. (10)) in reconstructing the decadal and multidecadal 
oscillations of the global surface temperature since 1850. 
The model has forecasting capability also at the decadal scale 
because the two curves calibrated using the independent periods 
1850-1950 and 1950-2011 are synchronous to each other also at 
the decadal scale and are synchronous with the temperature 
modulation revealed by the 4-year smooth curve: the statistical 
divergence between the harmonic model reconstruction and the 
data have a standard deviation of ø = 0.15 °C, which is due to the 
large and fast ENSO related oscillations, while the divergence with 
the grey 4-year smooth curve of the temperature has a standard 
deviation of ø = 0.05 °C, as Table 2 reports. 

Fig. 5 shows that the IPCC warming projection since 2000 (at a 
rate of 2.3+0.6 °C/century plus a vertical error of + 0.1 °C ) does 
not agree with the observed temperature pattern since about 
2005-2006. On the contrary, the empirical model we propose, 
Eq. (10), appears to reasonably forecast the observed trending of 
the global surface temperature since 2000, which appears to have 
been almost steady: the error bars are calculated by taking into 
account both the statistical error of the model (~ +0.1 °C) 
(because, at the moment, the harmonic model includes only the 
decadal and multidecadal scales and, evidently, it is not supposed 
to reconstruct the fast ENSO related oscillations) plus the pro- 
jected anthropogenic net warming with a linear rate within the 
interval 0.5-1.3 °C/century, as discussed above. According our 
model, by 2050 the climate may warm by about 0.1-0.5 °C, 
which is significant less than the average 1.2 +0.4 °C projected 
by the IPCC. If multisecular natural cycles (which according to 
several authors have significantly contributed to the observed 
1700-2000 warming and very likely will contribute to a cooling 
since the 21st century) are ignored, the temperature may warm 
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Fig. 5. [A] The monthly sampled global surface temperature from 1950 to 2050 (red); the proposed empirical model (Eq. (10)) made of the discussed 2 cycles (20 and 60 year) 
plus the quadratic trend until 2000 that is substituted with the corrected anthropogenic net projected warming as explained in the text (black); the IPCC (2007) projections 
(green). [B] The monthly sampled global surface temperature from 1950 to 2050 (red); a 4-year moving average estimates of the same (smooth wide gray curve); the proposed 
empirical model (Eq. (10)) made of the discussed 4 cycles (9.07, 10.44, 20 and 60 year) plus the quadratic trend until 2000 that is substituted with the anthropogenic net 
estimated contribution given by a linear trend with a rate within the interval 0.5-1.3 °C/century as discussed in the text (black and blue small curves); finally, by comparison the 
IPCC projected warming using the average GCM projection with a trend of 2.3 + 0.6 °C/century from 2000 to 2050. Note that the two harmonic model curves use the two decadal 
harmonics at 9.07-year and 10.44-year periods calibrated on the temperature data during two complementary time periods, 1850-1950 and 1950-2011 respectively. As evident 
in the figure, the decadal oscillations reconstructed by the two alternative models are very well synchronized between them and with the oscillations revealed in the grey 4-year 
smooth temperature grey curve. This validation result suggests that the astronomical harmonic model has forecast capability. The insert figure is reproduced in a full page figure 
in the supplement file. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 


by about 0.3-1.2 °C by 2100 contrary to the 1.0-3.6 °C warming 
projected by the IPCC (2007) according to its various emission 
scenarios. 

The divergence of the temperature data from the IPCC pro- 
jections and their persistent convergence with the astrono- 
mical harmonic model can be calculated by evaluating a time 
continuous discrepancy y?(t) (chi-squared) function defined as 


(Tem(t)—Mod(t))? 
(AMod(t))? 


Vib= (11) 


where Tem(t) is the 4-year smooth average temperature curve 
depicted in the figure, which highlights the decadal oscillation, 


Mod(t) is used first for indicating the IPCC GCM average projection 
curve and second for indicating the harmonic model average 
forecast curve as depicted in the figure, and AMod(t) is used to 
indicate the time dependent uncertainty first of the IPCC projec- 
tion and second of the harmonic model, respectively, which are 
depicted in the two shadow regions in Fig. 5. In the above 
equation the implicit error associated to the 4-year smooth 
average temperature curve is considered negligible (it has an 
order of magnitude of 0.01 °C) compared to the uncertainty of the 
models AMod(t), which has an order of magnitude of 0.1 °C and 
above, so we can ignore it in the denominator of Eq. (11). Values 
of y?(t) < 1 indicate a sufficient agreement between the data and 
the model at the particular time t, while values of y(t)>1 
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Fig. 6. The curves are produced with Eq. (11) and give a time dependent estimate of how well the astronomical harmonic model (crosses) and the IPCC projections (circles) 
forecast the temperature data (4-year smooth average temperature gray curve in Fig. 5) since 2000. The IPCC projections significantly diverge from the temperature data 
since 2004-2006. The astronomical harmonic model is shown to forecast the data quite well and it is quite stable in time. The variable y? needs to be less than 1 for 


statistical compatibility. 


indicate disagreement. Fig. 6 depicts Eq. (11) and clearly shows 
that the astronomical harmonic model forecast is quite accurate 
as the time progress since 2000. Indeed, the performance of our 
geometrical model is always superior than the IPCC projections. 
The IPCC (2007) projections significantly diverge from the data 
since 2004-2006. 


5. Discussion and conclusion 


The scientific method requires that a physical model fulfils two 
conditions: it has to reconstruct and predict (or forecast) physical 
observations. Herein, we have found that the GCMs used by the 
IPCC (2007) seriously fail to properly reconstruct even the large 
multidecadal oscillations found in the global surface temperature 
which have a climatic meaning. Consequently, the IPCC projec- 
tions for the 21st century cannot be trusted. On the contrary, the 
astronomical empirical harmonic model proposed in Scafetta 
(2010b, 2011b) has been shown to be capable of reconstructing 
and, more importantly, forecasting the decadal and multidecadal 
oscillations found in the global surface temperature with a 
sufficiently good accuracy. Figs. 5 and 6 show that in 1950 it 
could have been possible to accurately forecast the decadal and 
multidecadal oscillations observed in the climate since 1950, 
which includes a steady/cooling trend from 2000 to 2011. Four 
major cycles have been detected and used herein with period of 
9.1 year (which appears to be linked to a Soli/Lunar tidal cycle), 
and of 10-10.5, 20-21 and 60-61 year (which appears to be in 
phase with the gravitational cycles of Jupiter and Saturn that can 
also modulate the solar cycles at the equivalent time-scales). 
However, other astronomical cycles may be involved in the 
process. 

This result argues in favor of a celestial origin of the climate 
oscillations and whose mechanisms were not included in the 
climate models adopted by the IPCC (2007). The harmonic 
interpretation of climate change also appears more reasonable 
than recent attempts of reproducing with GCMs some limited 
climate pattern such as the observed slight cooling from 1998 to 
2008 by claiming that it is a red-noise-like internal fluctuation of 


the climatic system (Meehl et al., 2011) or by carefully playing 
with the very large uncertainty in the climate sensitivity to CO2 
changes and in the aerosol forcing (Kaufmann et al., 2011). In fact, 
a quasi 60-year cycle in the climate system has been observed for 
centuries and millennia in several independent records, as 
explained in Introduction. 

By not properly reconstructing the 20-year and 60-year 
natural cycles we found that the IPCC GCMs have seriously 
overestimated also the magnitude of the anthropogenic contribu- 
tion to the recent global warming. Indeed, other independent 
studies have found serious incompatibilities between the IPCC 
climate models and the actual observations and reached the same 
conclusion. For example, Douglass et al. (2007) showed that there 
is a large discrepancy between observed tropospheric tempera- 
ture trends and the IPCC climate model predictions from January 
1979 to December 2004: GCM ensemble mean simulations show 
that the increased CO, concentration should have produced an 
increase in the tropical warming trend with altitude, but balloon 
and satellite observations do not show any increase (Singer, 
2011). Spencer and Braswell (2011) have showed that there is a 
large discrepancy between the satellite observations and the 
behavior of the IPCC climate models on how the Earth loses 
energy as the surface temperature changes. Both studies imply 
that the modeled climate sensitivity to CO, is largely overesti- 
mated by the IPCC models. Our findings would be consistent with 
the above results too and would imply a climate sensitivity to COz 
doubling much lower than the IPCC’s proposal of 1.5-4.5 °C. 
Lindzen and Choi (2011) has argued for a climate sensitivity to 
a CO» doubling of 0.5-1.3 °C by using variations in Earth’s radiant 
energy balance as measured by satellites. We claim that the 
reason of the discrepancy between the model outcomes and the 
data is due to the fact that the current GCMs are missing major 
astronomical forcings related to the harmonies of the solar 
system and the physical/climatic mechanisms related to them 
(Scafetta, 2011b). 

Probably several solar and terrestrial mechanisms are involved 
in the process (Scafetta, 2009, 2010b, 2011b). It is reasonable that 
with their gravitational and magnetic fields, the planets can 
directly or indirectly modulate the solar activity, the heliosphere, 
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the solar wind and, ultimately, the terrestrial magnetosphere and 
ionosphere. In fact, planetary tides, as well as solar motion 
induced by planetary gravity may increase solar nuclear fusion 
rate (Grandpiere, 1996; Wolff and Patrone, 2010). Moreover, 
Charvatova et al. (1988), Komitov (2009), Mazzarella and 
Scafetta (2011) and Scafetta (in press) showed that the historical 
multisecular aurora record and some cosmogenic beryllium 
records presents a large quasi 60-year cycle which would suggest 
that the astronomical cycles regulated by Jupiter and Saturn are 
the primary indirect cause of the oscillations in the terrestrial 
ionosphere. Ogurtsov et al. (2002) have found that several multi- 
secular solar reconstructions do present a quasi 60-year cycle 
together with longer cycles. Loehle and Scafetta (2011) have 
argued that a quasi 60-year cycle may be present in the total 
solar irradiance (TSI) since 1850, although the exact reconstruc- 
tion of TSI is not currently possible. Indeed, TSI direct satellite 
measurements since 1978 have produced alternative composites 
such as the ACRIM (Willson and Mordvinov, 2003), which may 
present a pattern that would be compatible with a 60-year cycle. 
In fact, the ACRIM TSI satellite composite presents an increase 
from 1980 to 2002 and a decrease afterward. On the contrary, the 
PMOD TSI composite adopted by the IPCC Fröhlich (2006) does 
not present any patter resembling a 60-year modulations but a 
slightly decrease since 1980. However, the way how the PMOD 
science team has adjusted the TSI satellite records to obtain its 
composite may be erroneous (Scafetta and Willson, 2009; 
Scafetta, 2011a). 

Indeed, Scafetta (2011b) found that several mid-latitude aur- 
ora cycles (quasi 9.1, 10-10.5, 20-21 and 60-62 year cycles) 
correspond to the climate cycles herein detected. We believe that 
the oscillations found in the historical mid-latitude aurora record 
are quite important because reveal the existence of equivalent 
oscillations in the electric properties of the atmosphere, which 
can regulate the cloud system (Svensmark, 1998, 2007; Carslaw 
et al., 2002; Tinsley, 2008; Kirkby, 2007; Enghoff et al., 2011; 
Kirkby et al., 2011). In addition, the variations in solar activity also 
modulate the incoming cosmic ray flux that may lead to a cloud 
modulation. The letter too would modulate the terrestrial albedo 
with the same frequencies found in the solar system. As shown in 
Scafetta (2011b) just a 1-2% modulation of the albedo would be 
sufficient to reproduce the climatic signal at the surface, which is 
an amplitude compatible with the observations. Oscillations in 
the albedo would cause correspondent oscillations in the climate 
mostly through warming/cooling cycles induced in the ocean 
surface. For example, a 60-year modulation has been observed in 
the frequency of major hurricanes on the Atlantic ocean that has 
been associated to a 60-year cycle in the strength of the Atlantic 
Thermohaline Circulation (THC), which would also imply a similar 
oscillation in the Great Ocean Conveyor Belt (Gray and Klotzbach, 
2011). Moreover, herein we have found further evidences that the 
9.1-year cycle is linked to the Soli/Lunar tidal dynamics. Ulti- 
mately, the climate amplifies the effect of harmonic forcing 
through several internal feedback mechanisms, which ultimately 
tend to synchronize all climate oscillations with the solar- 
lunar-planetary astronomical oscillations through collective syn- 
chronization mechanisms (Pikovsky et al., 2001; Strogatz, 2009; 
Scafetta, 2010b). 

For the above reasons, it is very unlikely that the observed 
climatic oscillations are due only to an internal variability of the 
climate system that evolves independently of astronomical for- 
cings, as proposed by some authors (Latif et al., 2006; Meehl et al., 
2011). Indeed, the GCMs do not really reconstruct the actual 
observed oscillations at all temporal scales, nor they have ever 
been able to properly forecast them. It is evident that simply 
showing that a model is able to produce some kind of red-noise- 
like variability (as shown in the numerous GCM simulations 


depicted in the figures in the supplement file) is not enough to 
claim that the model has really modeled the observed dynamics 
of the climate. 

For the imminent future, the global climate may remain 
approximately steady until 2030-2040, as it has been observed 
from the 1940s to the 1970s because the 60-year climate cycle 
has entered into its cooling phase around 2000-2003, and this 
cooling will oppose the adverse effects of a realistic anthropo- 
genic global warming, as shown in Fig. 5. By using the same IPCC 
projected anthropogenic emissions our partial empirical harmo- 
nic model forecast a global warming by about 0.3-1.2 °C by 2100, 
contrary to the IPCC 1.0-3.6 °C projected warming. The climate 
may also further cool if additional natural secular and millennial 
cycles enter into their cooling phases. In fact, the current warm 
period may be part of a quasi millennial natural cycle, which is 
currently at its top as it was during the roman and medieval 
times, as can be deduced from climate records (Schulz and Paul, 
2002; Ljungqvist, 2010) and solar records covering the last 
millennia (Bard et al., 2000; Ogurtsov et al., 2002). Preliminary 
attempts to address this issue have been made by numerous 
authors as discussed in Introduction such as, for example, by 
Humlum et al. (2011), while a more detailed discussion is left to 
another paper. 


Appendix A. Supplementary data 


A large supplementary file with additional data associated 
with this article can be found in the online version at doi:10.1016/ 
j.jastp.2011.12.005. 
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